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' An extended Bose-Einstein condensate (BEC) in an optical lattice provides a kind of periodic 

\ dielectric and causes band gaps to occur in the spectrum of light propagating through it. We 

. examine the question whether these band gaps can modify the spontaneous emission rate of 

atoms excited from the BEC, and whether they can lead to a self-stabihzation of the BEC against 
spontaneous emission. We find that self-stabilization is not possible for BECs with a density in 
(~| ■ the order of 10^'' cm~'^. However, the corresponding non-Markovian behavior produces significant 

effects in the decay of excited atoms even for a homogeneous BEC interacting with a weak laser 
beam. These effects are caused by the occurrence of an avoided crossing in the photon (or rather 
polariton) spectrum. We also predict a new channel for spontaneous decay which arises from an 
^ I interference between periodically excited atoms and periodic photon modes. This new channel 

should also occur in ordinary periodic dielectrics. 
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Abstract 



1 Introduction 

It is well-known that the radiation properties of atoms can dramatically be manipulated by changing 
the environment where atoms emit photons. For micro-cavities it has been demonstrated and 
for periodic dielectric media predicted ||, || that a suppression of spontaneous emission (SE) can be 
achieved. In the case of a micro-cavity this happens because its geometry reduces the radiation-mode 
density, whereas in a periodic dielectric medium SE is suppressed due to the formation of photonic 
band gaps (PBG). 

The recent achievement of Bose-Einstein condensation in magnetic traps Q has provided a new 
state of matter where all atoms share a single macroscopic quantum state. Such a state of matter offers 
great opportunities to explore and test new phenomena related to macroscopic quantum coherence. 
Recently several authors have theoretically studied spontaneous emission in a trapped BEC. In this 
case the continuous center-of-mass momentum distribution leads to an increase of SE 10, |^. In 
addition, the stimulated emission can be increased by the Bose enhancement in a BEC In the 
case of two BECs, interference effects can be important [||. 

The present work is focused on the case of an extended BEC and was motivated by the following 
idea. If an extended BEC is placed in an optical lattice it will become periodic. Since such a BEC 
does provide a (quantum) dielectric it affects the properties of photons propagating through it and 
phenomena similar to PBGs do occur |1C)[. Since then the photon mode density around the resonance 
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frequency is reduced one can expect that SE is suppressed by non-Markovian effects. We thus are led 
to the foUowing question: can a BEC in an optical lattice stabilize itself against spontaneous emission? 

A large part of this paper is devoted to the answer of this question. However, we also have 
studied what will happen for a homogeneous extended BEC interacting with a weak running laser 
beam. Surprisingly, non-Markovian effects similar to that in a PBG do occur even in this non-periodic 
situation. This happens because in the presence of such a BEC photons and excited atoms do form 
superpositions called polaritons Q . The spectrum of these polaritons contains an avoided crossing 
which has a similar effect on the SE rate as a PBG. 

The paper is organized as follows. In Sec. ^ we will present the theoretical model on which our 
calculations are based. The general derivation of the SE rates in laser fields will be done in Sec. |^. 
The results for the case of a BEC in a traveling wave laser beam or in a ID optical lattice beam are 
discussed in Sees. ^ and || respectively, and are summarized in Sec. ^. The details of the calculations 
are given in two Appendices. 



2 The theoretical model 

We consider a BEC composed of two-level atoms which is coupled to the electromagnetic field. The 
interaction is described by using minimal coupling in rotating-wave approximation under neglection 
of the term quadratic in the electromagnetic field. The interaction Hamiltonian is then given by 

Hint = J dPkcPk'Q{k)a„{k)-^g{k')-^l{k + k')+YL.c. , (1) 

with Co-(fc) '■= oJi-csd-ea{k)[h/ {2{2n)^eoLLjk]^^^ for an electromagnetic mode with frequency ujk — c\k\ and 
polarization vector e^{k). The vector d denotes the atomic dipole moment and tOres is the resonance 
frequency. The Heisenberg equations of motion for the photon annihilation operators ao-(fc) and the 
field operators and for excited and ground-state atoms can be derived easily and are given by 




= '-^^g{k) + J d'k'Y,C{k'H{k')^e{k + k') (3) 
iM^ik) = huJka^{k) + Uk) f d''k'^l{k')^,{k + k') (4) 



We have neglected the interatomic interaction terms. 

To address the question of self-stabilization consider the following situation: the atoms in the 
internal ground-state have formed a BEC which is described by a macroscopically occupied coherent 
collective wavefunction They interact with a traveling wave or standing wave laser which is 

described by a coherent c- number field a'^^{k). Due to this interaction a part of the BEC is coherently 
excited. We denote the wavefunction for coherently excited atoms by 4'™'^. Since both the ground- 
state BEC and the laser beam are described by c-number fields it is easy to see from Eq. (||) that 
\['g°^ must be a c-number field, too. It is only through the spontaneous decay of these coherently 
excited atoms that q-number deviations from c-number solutions to Eqs. (H) to can appear. The 
corresponding SE rate determines the stability of the macroscopic solution. 

Let us start with the assumption that the BEC can indeed stabilize itself against SE. In that 
case a stationary macroscopic solution (5'™'^ ^g"*^, a^°'') of Eqs. (||) to (^) should exist. The problem 
then can be divided into two separate parts. We first search for the stationary macroscopic coherent 
solution which includes all interaction effects between atoms and photons beside SE. Having found 
this solution we can perform a stability analysis to analyze the quantum fluctuations (SE) around 
it. Spontaneous decay will make the coherent solution unstable and the corresponding quantum 
corrections will become important on a time scale comparable to the atomic lifetime (which is to be 
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calculated). For times shorter than this lifetime the deviations from the coherent solution will be 
small (i.e., there are only few non-condensed atoms and non-laser photons). 

Given a stationary macroscopic solution of the Heisenberg equations of motion the stability analysis 
can be performed by applying Bogoliubov's method. This is done by writing the quantum field 
operators in the form 

*,(fc) = exp[-ifit]{%°\k) + 6^gik)} (5) 

*e(fc) = exphl(Ai + WL)t]{^'™'^(fc)+(5*e(fc)} (6) 

a,(fc) = eM-i^Lt]{al°''ik)+Sa^ik)} (7) 

and retaining in Eqs. (^) to only terms linear in S'^i and 6aa, which describe the quantum 
fluctuations around the coherent solution. 

The resulting linearized equations of motions are given by 



iM*e(fc) = -nALS^e{k)+ J d^k'%°'^{k~k')Sa{k')C{k') (8) 

ihSifgik) = j Sk'-^l°^{k + k')5a){k')C{k') (9) 

iMd(fc) = h{c\k\-L0L)5a{k) + C{k) j Sk' \^^l°^\k + k')5-^l{k')+5-^e{k + k')-^f^*{k^^ 

Here ujl is the laser's frequency and := ojl — ^rca its detuning. The linearized equations are valid 
as long as the photon-atom quantum fluctuations remain small enough, i.e., there are only few non- 
condensed atoms and non-laser photons. This is certainly the the case for short times. Furthermore 
several other approximations have been made. First, it is not difficult to see that for a BEG of density 
lO^"* cm"'^ the chemical potential Ti^^ the kinetic energy Ti^k^/{2M) of an atom, and the laser's Rabi 
frequency Vl = a'^°''Ccro(^L)/?i are typically much smaller than the interaction energy C<T(^)^g°''(^) if 
I A: I is of the order of Wrcs/c We thus have neglected all terms in which these quantities do appear. 

In addition, we have introduced two specific polarization vectors for the electromagnetic field. A 
"non-coupled" polarization vector £Nc(fc)j which is perpendicular to the photon momentum k and the 
atomic dipole moment rf, and a coupled polarization vector ec{k)^ which is perpendicular to k and 
£*Nc(fc) (see Fig. Since the interaction is proportional to the scalar product of the polarization 
vector and d only electromagnetic modes with polarization ec{k) do interact with the atoms. We 
associate with these modes the quantum fluctuation operator 5a{k) :— dacr=c'{k)- It is easy to see 
that the scalar product ec{k) ■ d, which appears in the definition of (a-ik), is given by \d\ sin'd^, where 
is the angle between k and d. For notational convenience we have defined (^(k) :— Ca^cik)- 
It is possible to derive Eqs. (||) to (^ from an effective Hamiltonian for the quantum fluctuations, 

^^fiuct = Hpol + iJspont ■ (11) 

The first part, 

Hpol = h J d^k^^-ALS^lS-i!e + {c\k\-u;L)Sa''Sa^ 

conserves the number of photons plus excited atoms, 

Npoi = J d^fc{(5*^(5*e + Sa^Sa} . (13) 

The first integral in Hpoi describes the energy of free incoherent photons and atoms. The second inte- 
gral represents the excitation of atoms from the ground-state BEG and the reabsorption of incoherent 
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photons by the BEC. Its cigenmodes \q,r) — 7^1 ^|0) are generaUy superpositions of photons and ex- 
cited atoms, i.e., polaritons They are characterized by a continuous, momentum-like quantum 
number q and discrete quantum numbers r (see below) and can generally be written as 



q.r 



Sa{k) = 



d^k{£g-r{k)S'fl{k) + Ag,r{k)Sa\k)} 



r 

d^qY,^qAk)'Pq,r 



(14) 
(15) 

(16) 



The form of the expansion coefficients Sq'^rik), Aq'^rik) depends on the particular physical situation 
and is derived for a traveling and standing-wave laser in Appendix A. 2 and B.2, respectively. The 
second part of the effective Hamiltonian is given by 



spoilt 



^kd^k'C{k')-^l°'^{k + k') ^^Sa'' {k')S^l{k) + 5a{k')S^g{k)^ 



(17) 



It does not conserve A^poi and describes the spontaneous decay of coherently excited atoms. If this 
term would vanish the macroscopic coherent state would be stable against spontaneous decay. 



3 General derivation of SE rates 

The stability analysis essentially comprises to solve the time evolution of the polariton modes for 
relatively short times during which the occupation of the macroscopic coherent solution does not 
change very much. This will allow us to derive the initial SE rate of coherently excited atoms. We 
assume that initially all atoms and photons are in the state determined by the macroscopic coherent 
solution, or in other words, the quantized polariton field (photon-atom quantum fluctuations) is 
initially in the vacuum |0). This state then evolves under the action of the fluctuation Hamiltonian 
(11) into the time dependent state \ip{t)). 

To describe this time evolution we rewrite the polariton Hamiltonian ( p^ in the convenient form 

Hpol = f d^qJ2H^<l,r - ^L)Virrg,r , (18) 

r 

where u^^r — are the eigenfrequencies of Hpoi- Using Eqs. ( p^ ) and (^6[) one also can derive 

Hspont ^ Jd'kJ d^qY,{^'^9ik)'Pg-,rggAk) + S^lik)pl^gt^{k)] (19) 

with 

gg~r{k) / d'k'C(k')^l°'^{k + k')AgAk') ■ (20) 



To describe the evolution of the state \ip{t)) we make the following ansatz, which corresponds to 
the one-photon approximation, 

\m) « Ritm + Jd'kJ d\j2sqAk't)vy^i{k) \o) . (21) 

The Schrodinger equation ihlip) — -fffluctlV') then can be solved by using the Laplace transform R{s) = 
exp[^ts]R{t)dt and similarly for S^^sik,t). The resulting equations, 

ih{sR{s) ~ R{0)) = J d^k J d\J2^^Ak,s)ggAk) (22) 

ihsS^^k, s) = hiujq^r - ^LiSq^k, s) + R{s)g*^Ak) , (23) 



4 



have the solution 

The dependence on the particular physical situation is completely determined by the integral 

ih^ J ^ Zs- UJq^r ' ^ ' 

For notational convenience we have defined the complex variable 

Zs := IS + Ai . (26) 

The most important aspect of / is its complex analytical structure. This is because the inverse 
Laplace transform is defined by 
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f e + zoo 

m = 7h I e''R{s)ds , (27) 



where e is chosen so that the path of integration lies to the right of any branch cuts and poles of R{s). 
From Eq. ( |24| ) it becomes clear that the branch cuts of R{s) are those of / and that the poles of R{s) 
essentially depend on the form of /. Assuming that all poles of R{s) are simple poles we then find for 

m 

=y e*"'Res(i?(s),s,) + y"7r- / e*'R{s)ds , (28) 

where Si denote the poles of R{s) and Bj the branch cuts. Each pole corresponds to a fraction of the 
coherently excited atoms which decays (or increases exponentially) with a SE rate of 

7, = -2Re(sO- (29) 

The integration contours around the branch cuts corresponds to a fraction of coherently excited atoms 
with a non-exponential time evolution. 



With Eq. (|29[ ) we have found a general expression for the SE rates that can appear in the presence 
of a BEC. We now want to study the different physical situations of a BEC in a traveling or standing 
wave laser and to derive the corresponding values of 7^. To do so we have to find closed expressions 
for the polariton eigenfrequencies Ug-^r and the functions gq.rik) in order to derive I{s). These quan- 
tities in turn require the knowledge of both the polariton eigenmodes and the fields {'^g°^, O''^^^) 
comprising the macroscopic coherent solution. Since the calculations leading to a closed expression 
for / are quite involved we present them in the appendices. In the next two sections we analyse the 
results and give physical interpretations of the effects involved. 



4 Spontaneous emission rates for a BEC in a traveling wave 
laser 

In the case of a BEC interacting with a traveling wave laser, the polariton dispersion relation derived 
in Appendix |A.2| does contain an avoided crossing around the resonance frequency of the atoms (see 
Fig. y). Since nearly resonant photons provide the dominant contribution to SE, it is physically 
evident that this avoided crossing in the dispersion spectrum will produce an effect on the SE rate 
which is similar to what a PBG can do. 

This effect can be studied by analysing the closed expression for the renormalized value of I{s) 
which we have derived using the generalized Wigner-Weisskopf approximation presented in Appendix 



A.3. We find 



rRcn _ I W\ jRcn , ^e7vac j 47j/g 



5zs J biri 1 15z 



| + + (l-£i)3/Wh(./^)j 

(30) 
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In this expression 7vac := (?Wres/(37r?ieoc^) denotes the SE rate in free space and := Vpe the 
number of coherently excited atoms (V is the quantization volume). The frequency 

with pg being the density of atoms in the ground-state, determines the strength of the interaction 
between photons and excited atoms mediated by the BEC. Typically we have -Bintcraction = Ti^fgWrcs 
(see appendix |A.2| ). 

As already mentioned in Sec. ^ the time evolution of the macroscopic coherent solution depends 
essentially on the analytical structure of R{s) of Eq. (p4[). In general, R{s) has several poles and a 
branch cut originating from the term including the arcoth in This cut lies between Zg = and 

Zs — t^g- Another important property of is that it depends on Zg only through the ratio Zs/vg 
so that the magnitude of the SE modification depends on this ratio, too. In addition, a numerical 
evaluation of Eq. (^) shows that it is a slowly varying function of the order '^vacNe unless Zsjvg is 
close to zero. This has the following consequences. 

In free space z^ can be taken to be close to the pole of ^(s), i.e., of the order of jvncNe- Hence, 
the magnitude of non-Markovian effects is essentially determined by the ratio 

Zs_ _ iVe7vac _ 167r^ 

Vg^ Vg 3 PgA3 ■ 

This ratio is proportional to the total number of excited atoms divided by the number of condensed 
ground-state atoms per optical wavelength A. The request that this ratio should be small has important 
consequences when applied to BECs of a density in the order of 10^'' cm"^. In this case the BEC- 
induced effects can only be relevant if there are very few coherently excited atoms (in the order of 
one). However, for higher densities significant effects can occur also for a higher number of excited 
atoms. 

The request that is of order one also implies | | ^ VIl since otherwise the Rabi frequency 



would be large enough to excite many atoms. As discussed in Appendix A.l, the macroscopic coherent 
solution implies in this case for > the additional constraint A^ > Avg . In Fig. ^ the real part 
of the two dominating poles si,S2 of /(s) is shown as a function of A^ for the case A^e = 1 and 
— 2.57vac (corresponding to an atom density pg of 5 x 10^'' cm"'^). For A^ > a third pole 
appears with a very small negative decay rate (< 10~'^7vac)- The occurrence of negative decay rates 
is consistent within the range of validity of the linearized equations for the quantum fluctuations and 
may indicate the formation of an atom-photon bound state JlSt . Obviously the change in the SE rate 
can be quite large for small |Al|. According to Eq. ( p8| ) the fraction of atoms belonging to the poles 
can be easily calculated by determining the residue at the poles. It turns out that the pole whose 
real part asymptotically approaches 7vac always dominates and that the fraction of atoms belonging 
to other poles is significant only for small |Al|. The same holds for the fraction corresponding to the 
branch cut. 

If \ ^ 7vac holds the dominant pole si can be calculated by perturbation theory. Its real part 
(the decay rate) is then given by 



We see that the SE rate is altered by a factor of 1 — Ai/g/ {5Al)- It depends on the sign of Al whether 
SE is increased or decreased. We remark that the reason why SE depends on the detuning is that 
the coherently excited atoms are driven by the laser field and thus oscillate at the laser frequency lol 
instead of the resonance frequency ij^rcs (see Eq. (^). 

We shortly summarize the results that we have found for a BEC interacting with a running laser 
beam. For evolution times smaller than the atomic lifetime and for a weak laser beam which only 
excites a number = 0{1) of excited atoms (for BEC densities of 10^'* cm~'^), the spontaneous 
emission rate is significantly modified by non-Markovian effects. These effects result from an avoided 
crossing in the polariton spectrum caused by the extended homogeneous BEC in the internal ground 
state. 
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In the next section we will examine the corresponding results for a spatially periodic (lattice) BEC 
in a standing wave laser beam. 



5 Spontaneous emission of a BEC in a standing wave laser 



A BEC in an optical lattice formed by a standing wave laser beam becomes spatially periodic and 
thus provides a kind of periodic dielectric. The polariton spectrum wil l therefore contain band gaps. 
We have derived the corresponding dispersion relation in Appendix B_^. As is well known from PBGs 
a band gap around the resonance frequency will lead to non-Markovian effects in the spontaneous 
emission of an atom. The examination of these effects will answer the question whether a periodic 
BEC in an optical lattice can stabilize itself against SE. 

To determine the SE rate of a BEC in a ID optical lattice we have aga in to evaluate the integral 
/ of Eq. (H^). This task is quite involved and is presented in appendix B^. Our final analytical form 
of the renormalized integral J^''^" is given by the somewhat lengthy expression 



2711 I V^rcs 



vdv In 



vdv In 



{h{z,)-v){l-V- h(Zs)) + h{Zs) 
V^^" l'(«2 + l)l/4_(^2_^)l/4_ 



Zs foiZs) - V 



/i(w)arctanh 

ilvncNe Vg y/vdv 



+ 1)1/4 



h{v) 



/i(w)arctanh 



(«2 - 1)1/4 



+ 1)1/4 „(1„^) 



71" Zs Jo fo{Zs) - V 

/i('y)arctanh 



1/2 



{V^ + 1)1/4 



h{v) 



/i(w)arctanh 



(34) 



In this result we have introduced a couple of new notations. For notational convenience we have 
defined 



fo{zs) ■■= 



2ckj 



2cfci Zs 



Zs 2ckL 



(35) 
(36) 



as well as the abbreviation h{v) := \/ foizs) — fi{zs)/ {fo{zs) — v). We also introduced a cut-off 
A « rrieC^/h to regularize the integral (me is the electron's mass). Two important physical quantities 
are given by 



:= 



V 



(27r)3 

V 
(2^ 



coh \2 



coh 

e,2m— 1 e,2m+l 



(37) 
(38) 



where V denotes the quantization volume and the sum runs over the (real) momentum components 
^e°m of coherently excited atoms. Ng is simply the total number of excited atoms in the macroscopic 
coherent field, and iVg describes how these atoms are distributed in momentum space and is always 
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smaller than iVg. It is a measure for the degree of periodicity of the density of excited atoms, very 
roughly we have Vpe{z) ^ Ne + Ne cos{2zkL)- 

The influence of the BEC in a standing wave laser on the SE rate is determined by the frequencies 

•= T2 PaK^^I ~ 7^ (39) 

^3 •= T2 PaK^^I ~ . (40) 

In the polariton dispersion relation Vg produces a contribution similar to that of Vg in the case of a 
BEC in a traveling wave laser beam (avoided crossing). Vg produces a real PBG close to the resonance 
frequency due to the spatial periodicity of a BEC in a ID optical lattice. The two frequencies define 
the strength of the interaction mediated by the mean density pg and the periodic part pg of the 
ground-state BEC, 



Pa 
Pa 



C27]-)3 a,^ni) V*^) 



(2 



^E*^t™*™2™+2. (42) 



These densities play a similar role to what N^/V and Nf./V do for coherently excited atoms. For a 
mean density of pg ~ 10^'* cm^"^ we find i/g w 4 x 10^ Hz. The magnitude of Vg can vary between Vg 
for a very strong optical potential and if the laser beam is switched off. 

Although Eq. ( p4[ ) has a complicated structure it allows to analyze the main features of /f'-c" and 
hence of the time evolution of the macroscopic coherent solution in the presence of (small) quantum 
fluctuations. It is even possible to estimate the influence of the band gap with some simple arguments. 



5.1 General structure of the result 

A very important feature of the integral is that all parts of are proportional to Ne'^vac or 
■^e7vac- In addition, it becomes obvious that /^®" depends on the ground-state BEC and on the 
complex variable Zg essentially through the ratios Zsjvg and Zs/vg. The only exception to this is the 
first term containing Zg in Eq. (|35|), but this term is negligible compared to wres and does only serve 
to keep book on which side of the branch cut Zs is placed (see remarks below Eq. ([76|)). 

These facts can be exploited to estimate under which circumstances the influence of the BEC on 
the SE rate is significant. Since for a BEC in a traveling laser beam the contribution of the poles of 
R[s) usually dominates (see Sec. ^), we will focus on this part. The denominator of R{s) is of the 
form s — /^'^"(s), see Eq.(p^. Since a numerical analysis of Eq. ( p4| ) shows that /^''"(s) is of the order 
of its pre-factors Nejvac or Nejvac unless Zs/vg or Zsjvg are small, a pole Si must be of the order of 
these pre-factors (if the detuning Al is not very large). In analogy to the case studied in Sec. ^ one 
can again infer that the magnitude of the BEC-induced effects essentially depends on ratios of the 
form Ne^v&c/vg, for instance. As in Sec. ^ this allows the conclusion that for a BEC with a density 
in the order of 10^^ cm"'^ non-Markovian effects can only be relevant if there are very few coherently 
excited atoms (in the order of one). 

In the case of a BEC in an optical lattice this restriction has additional implications: the Rabi 
frequency O of the coherent standing wave laser has to be very small since otherwise too many atoms 
would be excited. The number of excited atoms is approximately given by Ne ~ (f^/Ai)^7Vg, where 
Ng denotes the total number of condensed ground-state atoms. Since Ng is very large the ratio fi/Ai 
must be very small in order to achieve TVe « 1. This, in turn, means that the optical potential 
(oc f2^/Ai) provided by the standing laser beam is very weak and thus Vg is much smaller than Vg. 
A small value of Vg simply means that the polaritonic band gap that is formed in the presence of a 
BEC will be small and therefore will not have a significant effect on the SE rate. 
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5.2 Interference channel for spontaneous emission in PBG 

Another observation deals with the dependence of I^"^ on the wavefunction of coherently excited 
atoms. It is known that in free space the shape of the spatial wavefunction of an excited atom does 
only have a tiny influence on its SE rate These small corrections are mainly due to the atomic 
kinetic energy which we have neglected in the Hamiltonian ( |lT|) for the quantum fluctuations. In 
this sense, one would expect that the SE rate in Eq. (jsj) does also not depend on the shape of the 
wavefunction for coherently excited atoms and is proportional to total number Ne of excited atoms 
in which this shape does not enter. However, Eq. ( |3^ ) does also include the terms proportional to the 
quantity Ng which clearly depend on the shape (for instance, vanishes for a spatially homogeneous 
wavefunction ^1°^). Principally, this contribution can be as large as that depending on N^. 

This new dependence on the shape of ^P^"'' is an additional effect of the BEC on the SE rate and 
no consequence of the polaritonic dispersion relation. To understand its origin it is useful to look 
at Eq. ( |6^ ) where iVg appears first. This contribution obviously does vanish if Ao{q,r) or A-i{q,r), 
which are momentum-components of the photon modes belonging to momentum q and q — 2fci, is 
zero. This is the case for photons interacting with a homogeneous BEC, for instance. The new effect 
therefore can be considered as arising from the interference between different momentum-components 
of the photon modes and the wavefunction of coherently excited atoms. 

We want to emphasize that this effect is not tied to the presence of a BEC. The only conditions for 
its existence are the periodicity of both the wavefunction 5*™^ for excited atoms and the eigenmodes 
for the photons. Since in an ordinary PBG material the photon eigenmodes are periodic, this new 
contribution to the SE rate can be present in ordinary PBG materials, too. In an ordinary periodic 
dielectric the new interference channel for SE even could produce large contributions since the period- 
icity of the dielectric is produced by, e.g., mechanical forces but not by the light that is used to excite 
the atoms. Only in the case of a BEC do the standing laser beams play a double role, excitation of 
atoms and production of a periodicity in the BEC, which results in a suppressed influence of both the 
polaritonic band gap and the interference channel on the SE rate. 

6 Conclusion 

We conclude this paper by summarizing the results that we have found. We have examined the self- 
stabilization of a BEC against SE by performing a stability analysis of a macroscopically occupied state 
for photons and two-level atoms, which describes a BEC that is coherently coupled to a laser beam. 
The presence of the ground-state BEC thereby leads to the formation of polaritons and introduces 
non-Markovian effects in the spontaneous decay of excited atoms. 

In the case of a BEC in a traveling-wave laser, the polariton spectrum displays an avoided crossing 
around the resonance frequency which causes similar changes in the SE rate as a PBG in periodic 
dielectrics. Its magnitude depends on the ratio Nf,/{pg\\) between the total number of excited atoms 
Nf. and the number of BEC-atoms inside a cube of the size of an optical wavelength Al. If this 
ratio is much larger than 1 the SE rate will essentially remain unchanged. Otherwise the change can 
be significant as the numerical examples shown in Fig. ^ demonstrate. The change of the SE rate 
displayed in Fig. || depends on the detuning of the laser because the coherently excited atoms are 
driven at the laser's frequency ujl- 

For a BEC in a ID optical lattice two new effects do appear. Being a kind of periodic dielectric 
the BEC then produces a real polaritonic band gap. The size of this band gap is determined by Vg. 
As in the case of a traveling wave laser, SE is only significantly altered if there are very few excited 
atoms. This in turn does imply that the optical lattice must be very weak and therefore produces 
only a small band gap which has only a very small influence on the SE rate. The second new effect in 
a periodic BEC is the appearance of a new channel for SE which arises from the interference between 
different momentum components of the excited-state wavefunction and the photon modes. Though 
its effect in a BEC is as small as that of the band gap it should also be present in the case of a PBG 
in an ordinary periodic dielectric where it can be large. 

It should be pointed out clearly where exactly the difference between an ordinary periodic dielectric 
and a BEC in an optical lattice comes into play. In an ordinary dielectric medium the periodicity is 
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produced by whatever forces determine the stabiUty of the medium. The excitation of an atom inside 
such a medium is done by a Ught beam, i.e., a completely different physical system. In the case of a 
BEC in an optical lattice, however, the excitation of the atoms and the potential that produces the 
periodicity of the BEC both are provided by the same device: the laser beams of the optical lattice. 
These lattice beams have to achieve two competing goals: to provide a strong periodic potential (to 
produce a large band gap) and to cause only a weak excitation (to have few excited atoms). As the 
achievement of both goals is impossible the periodicity of the BEC will only have a tiny influence on 
the SE rate and it will essentially cause the same effect as a homogeneous BEC in a running laser 
wave. 

This argument also provides the answer to the question whether self-stabilization of a BEC against 
SE is possible. Since the SE rate is only significantly changed if there are very few excited atoms, and 
since a large PBG does only form for strong laser beams, a self-stabilization is not possible for BECs 
with a density in the order of 10^^ cm~^. 

We finally remark that our results are not applicable to BECs confined in a micrometer-sized 
trap, a case discussed in the literature ^ 0]. Our work is concerned with BECs which are extended 
enough to allow the formation of polaritons. The necessary extension of the BEC can be estimated by 
considering the typical interaction energy for the formation of polaritons which is given by fi,JT7^uJ^ 
(see appendix |A.2| ). For a BEC with a density of 10^* cm~'^ this energy is in the order of ?i x 10^^ 
Hz. For the formation of polaritons a photon must therefore be inside the BEC longer than 10~^^ 
seconds. Since it travels at the speed of light the BEC must therefore be larger than about 3 mm. 

Acknowledgment: This work has been supported by the Australian Research Council and the 
Optik Zentrum Konstanz. 



A BEC in a running laser wave 

A.l Derivation of the macroscopic coherent solution 

We are interested in finding a particular solution (^'^°'\ a™'^(fc)) of macroscopically occupied 
fields to Eqs. (||)-(P which describes a BEC coherently coupled to a running laser wave. We thus 
make the ansatz 

(2^)3/2^5(fc)exp[-iAii] (43) 
exphic^it]J(fc - fci)5,,,,r!i[2(27r)3feoc^,Ji/2/(|d-|^^.^j (44) 

(2^)3/2^5(^ _ 4) exphi(/x + UL)t] (45) 

which corresponds to a homogeneous ground-state BEC of density pg, a laser beam with frequency 
wl, Rabi frequency il^ > 0, polarization ul, and wave- vector — k^ez {inside the BEC), and 
coherently excited atoms of density pe and of momentum TikL- Inserting these expressions into the 
Heisenberg equations of motions leads to a set of algebraical conditions which fix the chemical potential 
hp,^ the laser wavenumber k^, and the density of coherently excited atoms pe which we assume 
to be smaller than pg. If we neglect the kinetic energy the density of excited atoms is given by 
a/Pc = ^/Pg^L/ip + AL). The wavenumber kL is fixed by c/cl = ujl/2 + y/{ujL/2y - uj'^^^Vg/{p + Al), 
where Vg is defined in Eq. (^l|). Note that /cl generally is different from the free-space value ujl/c. 
For Al < (A^ > 0) the chemical potential is given hy p = — Al/2± {Al/2)^ + which implies 
p + Al > {p + Al < 0), respectively. Note that for > the expression for cfc^ implies the 
additional constraint p + Al > 4fg. 



a-^fc) = 



A. 2 Derivation of polariton eigenmodes 

Having found the macroscopic coherent solution we are in the position to to derive the polariton eigen- 
modes Because of the delta distribution appearing in the macroscopic solution (|4^) the Hamiltonian 
(12) reduces to a sum of two- level systems so that its eigenmodes are quite easy to find. They consist 
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of polaritons with momentum fiq and frequency spectrum 



^,-± ^^±W,, (46) 



where we have defined Wq := \J {Aq/2y + VgUj-ccs sin^ "dq and c|<7| — cjros- The coefficients of the 

polariton creation operator in Eq. ( p^ are given by 



/tjS + t/„CJres Sin^ 



so that we find for Eq. (^2 



The polariton spectrum lo±_^ clearly exhibits an avoided crossing around = of width 
^VgUj-rcs sin (see Fig. |2|). It also contains a small gap whose edge is reached in the limit ^ 
and oo Q. This gap will play no role for the SE rate since its edges are far away from the resonance 
frequency. This is physically reasonable since far away from resonance the two-level approximation 
for the atomic internal structure ceases to be valid. 

A. 3 Calculation of I{s) 

With the help of Eq. ( |48| ) the integral /(s) of Eq. ( [25| ) can be written in the form 



f27r)3 I " '" ( . . i/g litres sin^ iJ- 

where V denotes the quantization volume. This integral agrees with the one found in absence of a 
BEC (which describes SE in free space) by setting Vg = A^ = 0. We denote this free space integral 
by Iq :— lii'g = A/^ = 0). Both integrals are linearly divergent and can be treated in the way pointed 
out by Bethe (see, e.g., Ref. |l^), i.e., we renormalize the integrals by subtracting the free-electron 
contribution, 

I^-^-.^I-Y;^ fd^k'^ (50) 
(27r)3i J ul ^ ' 

and 7^''" = I^™{yg = Al = 0). These renormalized integrals are only logarithmically divergent. 

At this point it is customary in the calculation of the free-space SE to perform the Wigner- 
Weisskopf approximation by neglecting the dependence of I^"^{s) on s. In the presence of a band 
gap this is inappropriate due to the strong variation of the mode density around the gap ||^, ^. 
Nevertheless, we can perform a generalized Wigner-Weisskopf approximation in the following way. 
We expect that the typical timescale on which SE happens is much larger than the optical cycle 
timescale l/wres- From the definition of the inverse Laplace transform ( p7| ) it is clear that the variable 
s plays more or less the role of a Fourier-transformed time. We thus expect that only values of s much 
smaller than Wrcs contribute significantly to the SE. This implies that we can neglect (the imaginary 
part of) s wherever it appears together with Wios ot ul- Thus, we are allowed to set s — ilul ~ —ii^L 
in the denominator of I{s) while retaining the term depending on s — iAi,. 

In the case of /(^°" this procedure immediately reproduces the Wigner-Weisskopf result 7^''" w 
^e{(7vac/2) -HiA2±y^}, where 7vac and are defined in Sec. |. To fix /^l^Zh follow the theory of 
Bethe (see, e.g., [|2|) and introduce a cut-off frequency of rrieC^ /Ti in 7^°^^, where is the electron's 
mass. Calculating the principal value of the integral then leads to A^a^i, ~ 27vac- In contrast to free 
space the SE rate in a BEC depends on A^a^^, since such a radiative frequency correction shifts the 
center of the avoided crossing (or of a band gap |^). 

It remains to calculate a renormalized expression of the integral 7^^^" in the presence of a BEC. 
Fortunately, this task reduces to integrals proportional to 7;^''" and a couple of convergent integrals 
and leads to Eq. (pG) 



11 



B BEC in a standing wave laser 



B.l Derivation of the macroscopic coherent solution 

Since in a running laser wave the BEC density is not periodic no PBGs can be formed. It is therefore 
of interest to study a BEC interacting with a standing laser wave so that the formation of photonic, 
or rather polaritonic band gaps is possible. The coherent laser field describing a standing wave 
is given by 

a^'^lfc) = al^'^S^^^.iSik-kL) + 6ik + kL)}e^p[-iLULt] . (51) 

We assume that the amplitude af^ is real and that the polarization (To of the laser beam is parallel 
to the dipole moment d of the atoms. 

Since the laser field provides a periodic potential for the atoms it is reasonable to assume that the 
macroscopic atomic fields are periodic, too (at least for the ground-state of the system). One also can 
make the ansatz that 'i'g°^ has period 2kL so that the coherent solutions can be written as 

M.l°\k) = e-*^*^5(fc-2nfci)vl/-h^ (52) 

n 

vf^^fc) = e-'('^+--)*^5(^-(2n+l)fcL)vI/-';,+i. (53) 

n 

Inserting this into Eqs. @ to (^) leads to the matrix equations 

iu:L-ckL)n = ^^i^^^coh^vl/coh^^^ (54) 

n 

( fiP\ 

Ai+A.-(2n+l)2-^U'™h^i = + *™2„+2} (55) 



(^/i-(2«)'^j*™2^. = m*™2Vi + *™2';+i}- (56) 

We have assumed that Cct(A;) is real and does not depend on the sign of k and introduced the real 
Rabi frequency :— a'^{'^C^„^^(ki}j jh. For consistency with the assumption that af^ is a real quantity 
the coefficients and must be real, too. 

The system (p4[)-(p6|) of algebraic equations can easily be solved numerically. To do so we assume 
that the Rabi frequency of the laser beam is a given quantity. For a given value of kj^ the two equations 



(55) and (£^) then just describe the well-known problem of a two- level atom moving in a standing 
laser wave. This is a simple system of linear equations and can be solved in a standard manner. The 
solution then can be inserted into Eq. (|5^ ) which then, because loj^ and a™^ are fixed, determines the 
value of fc^. We then have reinserted the new value for kj^ into the system (|4|)-(^6|) and iterated the 
procedure until fc^ did not change significantly anymore. 



B.2 Derivation of polariton eigenmodes 

In this case the periodic structure of the macroscopic solution ( p^ ) leads to a more complicated 
structure of the eigenmodes of iJpoi than in the case of a traveling laser wave. To find these modes 
we make in Eq. ( [1^ ) the ansatz E^^r(k) — X^mez '')'^(^ ^ 9^ 2mki^) and correspondingly for 

^g,r(fc)- Now q denotes the quasi-momentum of the polariton. For a single standing laser wave along 
the z-axis is confined to [— A;i,/ci] whereas qx^qy represent the real momentum of the polariton 
perpendicular to the laser beam. The index r is a collective notation for discrete quantum numbers 
which include an internal quantum number taking two values (since two quantum fields Sa and S'^e 
are involved) and the band index. This ansatz leads to 

9qAk) := J2 r)Ciq + 2m^L)*™h(^ ^ 2mkL) . (57) 
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and results in the eigenvalue equations 

E{q,r)£,,Mr) = -nAi£,„(g, r) + ^ ^^'^^.Cl? + 2(m ~ n)fci)An-„(g, r) (58) 

n 

E{q,r)Am{q,r) = h{c\q + 2mkL\ ~ u;L)Arniq.r) +J2%tnC{q + '2mkL)£rn+niq,r) . (59) 



These equations can be substantially simplified by noting that the frequency difference h{c\q+2mkL \ — 
lol) is huge compared to all other energy scales involved unless m — 0,±1 and qis close to ifc/,. We 
thus can approximate the photon-part of all modes with |to| > 1 as free photons and need only to 
retain the coefHcients Ao{q, r) and A-i{q, r) for q^ G [0, ki,] and the coefficients Ao{q, r) and Ai{q, r) 
for qz £ [— fcz,, 0], respectively. We will focus here on the case qz E [0, /cl] since the second case can be 
treated analogously. 

To solve the resulting equations we introduce the two quantities Fo :— J2n'^g°2n^n{q,r) and 

F^i :— ^n'^g°2n^n-i{q,r) and make the approximation ({q) w C{q ~ 2fcL) ~ C(^l) so that the 
problem is reduced to the simple matrix eigenvalue equation 



{E{q,r)+hAL) 



Here we have introduced 











Ao 




V -4-1 ) 







C(fcL) 



C(fcL)(27r)3pg ^(^^)(2;r)3p, \ 
C(fcL)(27r)3p3 C(fcL)(27r)3p<, 







hAo 








\ 


( ^" \ 








Ao 


J 


\A-, J 



(60) 



Ao c\q] - LJrcs 

A_i := c\q- 2kL\ - c^rc 



(61) 
(62) 



and the quantities Pg and Pg which are defined in Eqs. ( p] ) and (p^), respectively. 

The problem of finding the eigenvalues and eigenvectors of a 4 x 4 matrix is a basic one. The 
eigenvalues lo^,. := {E{q, r) + hA]^)/h fulfill the relation P^hi'-Oq.r) = 0, where 

Pch{z) [z^ - zA^i - VgUj,.cs){z'^ - zAo - DgUJi-cs) - i>g(^ycs (63) 

is the characteristic polynomial of the matrix and the frequencies Dg and Dg are defined in Eqs. (|39| ) 
and @. 

Due to the periodicity of the ground-state BEC function of the quasi-momentum q exhibits 

the phenomenon of band gaps , see Fig. ^. Though closed expressions for the eigenvalues ujq-^r do 
exist they are rather cumbersome and not of much use for our problem. 

We instead will use a theorem on the eigenvalues to derive the physical quantities of interest. To 
do so we assume that we already know the eigenvalues. For a given eigenvalue oj^^r it is easy to solve 
for the eigenvectors. For the relevant components we find 



-4o(g,r) = 



A-iiq,r) 



^2l>2tj2^_, + (^ujl^ + t'gW,.os)(2w|^ - Wg-,r(Ao -I- A_i) - 2t'gW,.os) 



(64) 



(65) 



^2D'^LJ^^^ + {uji^ + l^gUJrcs)i2ujl^ - a;^,r(Ao -I- A_i) - 2l^gUJrcs) 

The normalization has been done by requiring the particle number (|l3|) to be pseudo- normalized, i.e., 

Y,{\£rn{q,r)\' + \AUq,r)\'}^l. (66) 

m 

The analogous equation for q^ G [— fci,0] is obtained if A-i{q,r) is replaced by Ai{q,r) and A_i by 
Ai := c\q + 2kL\ - Wres- 
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B.3 Calculation of /(s) 



To calculate the mtegral I{s) of Eq. (25) we first need to calculate the integral over k in the nominator 



which is an easy task. With Eq. ( |S7| ) we find for € [0, fc^] 

d3fc|g-,(fc)|2 ^ N,{Al{q,r)e{q}+AlMr)e{q- 24)} + 2N,A-i{q,r)Aoiq,r)C{q)C{q~ 2fci 



(67) 

if r is in the lowest two energy bands. The numbers and Ne are defined in Eqs. ( p7| ) and (p^), 
respectively. For q^ e [— fci,0] and the two lowest energy bands one has to replace A-i{q,r) by 
Ai{q,r) and '2kL) by (^{q + 2/0^). For all higher bands, which according to our approximation 

just describe free photons, the integral ( |67|) has the simple value NeC,^{q). 

From Eq. (^) it becomes clear that to calculate / of Eq. ( p5| ) we have to find closed expressions 
for terms like A-i{q, r)Ao{q, r)/{zs — Wg*,r)- For the lowest two energy bands, the sum over r now 
runs over the four eigenvectors of the matrix (|60|). It would be extremely tedious if not practically 
impossible to derive these sums by simply inserting the complicated closed expressions for u!^,r into 
them. Instead, we start with the observation that the polynomial appearing in the denominator of 
Eqs. (^4|) and ( |65| ) can be written as 

where P^i^{z) denotes the derivative of the characteristic polynomial (|63|). This enables us to write 
the sum under consideration in the form 

^ A-i{q,r)Ao{q,r) ^ ^ uj^,r 



{Zs - UJq;r) P^y^{uJ^^r){Zs - UJq^r) 



This can be further simplified by noting that the characteristic polynomial can also be written in the 
form Pch{z) = Yiri^ ~ ^q.r) SO that we have f'ch('^g,r) = Y\r>^r^'^q,r ~ '^q,r')- Using this expression it 
is straightforward if still tedious to find 

■A-i{q,r)Ao{q,r) ^ Zs^gUj^cs 

^ {Zs-U!q^r) Pch{Zs) 

We thus have been able to calculate this sum without explicit knowledge of the eigenvalues of Eq. (j6^) . 
The remaining sums which result from the insertion of Eq. (^7|) into Eq. (|2^) can be treated in a similar 
way and are given by 

^ {Zs-UJq,r) PchiZs) ^ ' 

^-l(g,0 ^ ^ Z^ - Z.Aq- UgUres ^^2) 

for the lowest two energy bands and qz € [0, k^]. Again the corresponding expressions for q^ G [— fci, 0] 
are obtained by replacing A-i{q, s) by Ai{q, s) and A_i by Ai. 

Taking everything together the use of Eq. (|6^) and the three sums (^ , ( [n] ) , and (|7|) allows us 
to bring the integral ( |25| ) into the form 

/ = 



— 2 / dqz dq^dqy 



NeC{q)'^[zl - ZsAl - VgUJ.cs] + NeC,{q + 2kLf[zl - ZsAp - DgLOres] + 2Nei>gL0resCiq)Ciq + ^kL) ^ 

ih^ 7y_ zl - zf (Ao + Ai) + z^AqAi - 2j>gW,es) + 2;.i^gWrcs(Ao + Ai) + dIlo^^^ - vjuj^^^ 

~s I ^:^^NeCiq)^[Zs - ZsA_i - DgUl^cB] + NeCiq- '^^l^Iz^ - Z^Aq - i^gWres] + 2N eVgi^rcs(^{q)C{q - S^l) 



iT? Jv+ " zl- z3(Ao + A_i) + z2(AoA_i - 2DgLj,^,) + z,Dgu;,^,{Ao + A_i) 

(73) 
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lo = -pi I dqxdqydqz- ^| , — (74) 



where the areas of mtegration V± are given by qx,qv G (—00,00) and qz € [0, /c^] for V+ as weh 
as qz e [— fcijO] for V-. The first integral in Eq. ( |73| ) represents the contribution from the higher 
polaritonic energy bands where the polaritons can be considered as free photons. Apart from the 
restriction \qz\ > 2fci, which essentially means that the photons are far off-resonant, it has the same 
form as the integral 

—2 / dq.j;dqydqz 

i?i J-00 - c\q\ 

which appears in the calculation of the free-space SE rate. The second contribution arises from the two 
lowest energy bands for negative quasi-momentum > qz > —ki,, and the third integral represents 
the corresponding contribution for positive quasi momentum g^. It is not hard to see that Eq. (73) 
reduces to the free integral (^J) in absence of a ground-state BEC, i.e., for = t>g = and to Eq. (4^) 
if the BEC is homogeneous (for I'g = but D „ ^ 0) . 

We now return to the evaluation of Eq. ([73|). Shifting the integration variable qz and exploiting 
the symmetries C{~k) = C{k) and Ai{~q) — A^i{q) allows us to combine the last two integrals into a 
more convenient form. To simplify the process of renormalization it is also advantageous to calculate 
I — Iq instead of / alone since in this difference the divergence appearing in Iq is canceled. We then 
find 



—2- / dqz / dq^dqy — 

/"""^J ^ 7VeC((?)2[z2-Z,A_i-PgWres]+2iVeZ>gW,esC(9)C('f+2fcL) 



I-h = I dqz I dq^dqy- 1^^' ^ + (75) 



• 4-2 / / dqxdqy . „ _ _ 

ih Jo J [z; - ZsAo - VgUJr-cs] [zi - ZsA-i - l^gUJres] " l^g^^rcs 

To calculate these expressions we have to make one further approximation by neglecting the angular 
dependence in C(<z) = sin((?)ijji.cs|rf| [?i/(2(27r)"^eoc|(?|]^^^, where 6 is the angle between q and the atomic 
dipole moment d. Replacing sin(0) by •\/2/3 for all values of g leads to the correct result in absence of 
a BEC and should produce qualitatively correct results for the case under consideration. Doing this 
approximation in the case of a homogeneous BEC just amounts in replacing a factor of 4/5 by 2/3 in 
the modifications of the SE rate (p3|). We remark that neglecting the dependence of C, on sin(6') does 
only symmetrize the integrand in the [q^^qy) plane. This does not correspond to an isotropic band 
model because the asymmetry between qz and (g^r, qy) still persists. 

The calculation of the first integral of Eq. ([75|), which roughly corresponds to the contribution of 



— /o, is quickly done and results in 

—2- / dqz / dqxdqy — = . <^ 2 + 2 In + i7rsgn(Im(z5)) \ . (76) 

h Jo J -00 Zs - C\q\+ Wrcs 27ri |^ V^ics/ J 

Here A is a cut-off which usually is taken to be A = MeC^/?i, where Me is the electron's mass (see, 
e.g., Ref. 

The dependence on the sign of Im(zs) originally comes from a logarithm of the form ln(— (cjics + 
Zs)/A). This expression can be reduced to the one presented in Eq. ( |7^ ) by doing a generalized 
Wigner-Weisskopf approximation as it was introduced above in the case of a running laser wave. 
Though Im(zs) is also much smaller than any other quantity in the above logarithm, it determines the 
sign of the imaginary part of the logarithm's argument. Since the logarithm has a branch cut along 
the negative real axis, this sign determines on which side of the cut we are. 

It is also worth remarking that Eq. ( [76|) is logarithmically divergent with A although we did not 
subtract the free-electron part, a step which in the free-space calculation is done to remove a linearly 
divergent contribution (see, e.g., Ref. [|l^). This is because the integration over qz does not extend to 
infinity, thus reducing the degree of divergence by one. 

As has been already mentioned, Eq. (^6|) very roughly corresponds to the negative of the free- 
space integral /q. As a consequence, its contribution will be mostly canceled after the renormalization 
of /. This renormalization is easily done by noting that / — /q = /Ren _ /Ren ^j-^g^^ /Ron _^ 
{I — Iq) + /(^°", where the superscript "Ren" denotes the renormalized integrals and the renormalized 
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free-space integral is approximately given by 



rRcn ,7vacA^e ( A 



^Hcn ^ ^^v^ ^_ J ^ ,^3gn(i„i(^,))| . (77) 
This allows us to derive from Eq. ( |75| ) the expression 

,Rc„ _ ^e7vac f ^ + \n ( — ) } + (78) 



27ri V t^i 



Ain'^uJr-cs Jq ^ J ^ \q\ {[z^ - ZsAq - iygUj^cs][z^ - ZsA^i ~ VgUj,.cs] - v^uj^cs} 

' dqz / dqxdqy 



To reduce Eq. ( [78|) it is convenient to introduce the scaled variables of integration u :~ qz/{2kL) and 

i) 

Ren 



it is convenient to introduce the scaled variables of integrat: 
V :— {ql H- qy)/{4:kj^) and the abbreviations (|35| ) and ( [3^ ) for the evaluation. Eq. ([78D then becomes 



""^^^ J2 + ln(Au+ (79) 



7vac^e r {ckL/u;,,s)[fo{Zs) - ^ {u - If + v] 



2ni Jo Jo V^FT^{[fo{zs) - VI?T^][/o(z,) - ^{u - 1)2 +v]- /i(z,)} 



47ri Jo Jo [foizs) - V^^^][fo{zs) - V(« - 1)' + - hi^s) 



This expression can be further reduced by switching to the integration variable v' — {u— 1)^ + v 
and exchanging the sequence of integration so that one first integrates over u. This allows us to reduce 
the integral / to a number of one-dimensional integrals and leads us to our final analytical result (^). 
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Figure 1: The polarization vectors can be chosen in a way that only one of them, ec{k), is not 
orthogonal to d. In this case the three vectors d, k, and £c{k) are in the same plane. The second 
polarization vector £Nc(fe) is perpendicular to this plane. 




Figure 2: A homogeneous BEC induces an avoided crossing in the polariton spectrum. Far away 
from the avoided crossing the polaritons describe excited atoms or photons. Thus, if one focuses on 
the photons, the avoided crossing provides an effective band gap. 
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Figure 3: Spontaneous emission rate of a partially excited BEC in a running laser wave for detuning 
Ai > Ai/g (solid lines) and < (dashed lines). The atoms break up into different fractions with 
different decay rates. The dominating fraction is the one whose decay rate asymptotically approaches 

"Tvac* 
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